set term postscript eps
set output "q7.ps"
a=1
b=1
c=6
d=0
f(x) = a*exp(-b*(x-c)**2) + d
fit f(x) "q7.dat" u ($0+1):1 via a,b,c,d
plot "q7.dat" u ($0+1):1 w points, f(x)

